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Abstract 

We have applied our empirical-PSF-based photometric techniques on a large number of calibration-related 
WFC3/UVIS UV-B exposures of the core of cj Cen, and found a well-defined split in the bright part of the white-dwarf 
cooling sequence (WDCS). The redder sequence is more populated by a factor of ~ 2. We can explain the separation of 
the two sequences and their number ratio in terms of the He-normal and He-rich subpopulations that had been previously 
identified along the cluster main sequence. The blue WDCS is populated by the evolved stars of the He-normal compo- 
nent (~ O.55AI0 CO-core DA objects) while the red WDCS hosts the end-products of the He-rich population (~ 0.46 Al© 
objects, ~ 10 % CO-core and ~ 90 % He-core WDs). The He-core WDs correspond to He-rich stars that missed the 
central He-ignition, and we estimate their fraction by analyzing the population ratios along the cluster horizontal branch. 



Subject headings: globular clusters: individual (NGC 5139) — Hertzspmng-Russell and C-M diagrams 
white dwarfs — techniques: photometric 



stars: Population II 



1. Introduction 

The massive globular cluster (GC) w Cen has been studied 
since the '60s because of its peculiar stellar populations (see 
van Leeuwen, Hughes & Piotto 2002 for a review). The dis- 
covery of its multiple stellar generations (Bedin et al. 2004, 
Bellini et al. 2010) has further revitalized detailed photomet- 
ric and spectroscopic investigations of this object. The white- 
dwarf (WD) cooling sequence (CS) is the least-studied part of 
the color-magnitude diagram (CMD) of this cluster, for the ob- 
vious reason that WDs are faint and are usually concentrated 
in the most-crowded central regions. 

Ortolani & Rosino (1987) were the first to detect a dozen 
WDs in (L) Cen, but it is only with the advent of HST that it 
became possible to analyze the WDCS in more detail (Mon- 
elli et al. 2005; Calamida et al. 2008). In particular, Calamida 
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et al. showed that the WDCS spread in color is larger than 
expected from photometric errors alone. We will show here 
that the WDCS of <jj Cen actually splits into two distinct se- 
quences, which can be explained as separate WD populations 
of different mass and chemical stratification, corresponding to 
the evolved components of the separate subpopulations with 
differing helium abundances, as revealed by the photometric 
analysis of main-sequence stars (e.g. King et al. 2012 and ref- 
erences therein). 

2. Data sets and reduction 

This analysis is based on HST's WFC3/UVIS observations 
taken between 2009 and 2011 under WFC3 calibration pro- 
grams. Specifically: we have 27 x 900 s exposures in F225W, 
31 X 800 s in F275W, 37 x 350 s in F336W, and 34 x 350 s in 
F438W. The field of view (FoV) is centered on the cluster's 
center, and covers the inner ~ 150". 

All WFC3/UVIS _flt images for a given filter were an- 
alyzed simultaneously to generate an astrometric and photo- 
metric catalog of stars in the field, using an evolution of the 
software tools described in Anderson et al. (2008). Briefly, we 
used information from all the available exposures for a given 
filter to find significant detections across the FoV. These de- 
tections were then fit independently in each exposure by em- 
ploying a fixed array of 7 x 8 spatially-varying static hbrary 
PSF models, plus a 5 x 5 array of perturbation PSFs tailored 
to fit the star profiles in each exposure in order to account for 
focus/breathing variations. Positions were corrected for geo- 
metric distortion according to the routines in Bellini, Ander- 
son & Bedin (2011). The photometry was calibrated into the 
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Fig. 1 . — The double WDCS of lo Cen: from left to right we show mF438w versus mF225w -'nF438w, 'Kf^vsw -"2f438w. and raFsaew -"'F438W. respectively. 




Fig. 2. — (a) Blow-up of the WDCS in the raF275w versus raF275w - "'F438W CMD with the four points used in the coordinate-system transformation 
of the next panel, (b) Rectification of the WDCS region that was shown in gray in panel (a). Here we define as bWD (in blue) those stars having 
"abscissa"<0, rWD stars the others (in red), (c) The FoV is divided into three equally-populated radial intervals. For each of them we fit a dual 
Gaussian to the "abscissa" distribution, and used the area of the two Gaussians to estimate the number of bWD and rWD stars in each radial interval 
(radial extension and areas are quoted in each panel), (d) Radial distribution of the ratio of rWD with respect to the total WD population. 



Vega-mag flight system following Bedin et al. (2005), using 
the WFC3 zero points provided by the STScfl 

3. The Double WD CS 

Figure 1 shows the WDCS of o) Cen in three different 
CMDs, keeping the »zp438w magnitude on the vertical axis, and 
using a variety of color combinations. WDs are shown with 
heavier dots, for clarity. We identified over 2,000 WDs across 
the FoV. The WDCS can be followed for over 6 magnitudes in 
«^F438w, starting from just below the termination point of the 
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extreme horizontal branch (HB). Two distinct and well defined 
WD sequences are visible in the bright part of the CS, before 
becoming indistinguishable below otf438w ~ 25. A few stars, 
about 1 magnitude redder than the WDCS, define an additional, 
scarcely populated sequence. 

Stars on the main sequence (MS) of oj Cen are well known 
to be separated into three main components: (1) a He-rich 
bMS, (2) a He-normal, [Fe/H]-poor rMS, and (3) a [Fe/H]- 
rich MS-a (Anderson 1997; Bedin et al. 2004; Piotto et al. 
2005; Bellini et al. 2010). The bMS has been found to be 
more centrally concentrated than the rMS (Sollima et al. 2007; 
Bellini et al. 2009), a fossil signature of the formation process 
of the two populations. Indeed, the present-day relaxation time 
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(fieiax = 12.3 Gyr, Harris 1996, 2010 edition) suggests the two 
populations have not yet had enough time to mix. 

In order to look for a connection between the split WDCS 
and the split MS, we measured the radial distribution of the 
two WD components as follows. We started by selecting stars 
in the mpns^ - 'Wf438w CMD (Figure 2a). To optimally sepa- 
rate the populations, we used a procedure developed by Milone 
et al. (2009). We first selected by hand two points (Pib,P2b) on 
the blue CS, and two points (Pir,P2r) on the red one. We then 
linearly transformed the CMD into a reference frame in which 
the coordinates of Pib, P2b, Pir, and Pa,, map to (-1,1), (-1,0), 
(1,1) and (1,0), respectively (Figure 2b). Here the blue com- 
ponent (the bWD) maps to negative "abscissa" values, and the 
red component (the rWD) to positive "abscissa" values (color 
coded in blue and red, respectively, in panel (b)). To examine 
the radial distribution we followed the methods by Bellini et 
al. (2013). We divided the FoV into three equally-populated 
radial intervals, and fit in each of them a Gaussian to each of 
the color-like distribution of bWD and rWD stars (panels (c)). 
(A dual-Gaussian fit is more robust than simply adding up stars 
in histogram bins.) The radial coverage of interest and the 
areas of each of the two Gaussians are given in the top-left 
corner of each of these panels. Panel (d) shows the quantity 
PrWD = A^rWD / (MwD +M)Wd) as a function of distance from the 
cluster center Binomial errors are ~ 3.7 %. There is marginal 
evidence for a gradient, with rWD stars more concentrated that 
bWD ones. Moreover, rWD stars (326) are ~ 64 % of the total 
bWD-t-rWD stars (507). We note that Bellini et al. (2010) also 
did not find any strong evidence of a radial gradient in the ratio 
of the two main MSs in the inner 2 arcmin of w Cen (see their 
Fig. 18). 

The top and middle panels of Fig. 3 show the WDCS CMD 
in a wide variety of color systems, with the same stars color- 
coded the same ways in all the panels. The split is present in 
all of them. It is interesting to note, however, that even though 
all the other populations (MS, sub-giant branch, etc) show dra- 
matic splits in two-color diagrams (2CDs), there is no split in 
the 2CDs plotted in the bottom panels of Fig. 3. This implies 
that the split is likely due to a discrete mass/radius difference 
between the populations, rather than an atmospheric effect. 

4. The nature of the two WDCSs 

We have adopted the distance modulus derived from the 
cluster eclipsing binary OGLEGC 17 (Thompson et al. 2001), 
i.e., {m - M)o - 13.66 + 0.12, and a mean value of the red- 
dening E{B - y) = 0.12 as in King et al. (2012). To calcu- 
late the extinction in our filters we have applied the extinc- 
tion law of Cardelli et al. (1989, with Ry = 3.1) to our WD 
spectra used for the calculation of the bolometric corrections 
to the UVIS filters i The resulting extinctions in the F275W, 
F438W and F336W filters are: 0.75, 0.50 and 0.62 mag, re- 
spectively. These values were applied to our WD tracks - 
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Fig. 3. — Top and middle panels: multicolor analysis for the WDCS. 
Bottom panels: two-color diagrams. Red and blue WDs are those 
defined in Figure 2b. 



transformed to the UVIS filters - from BaSTI CO-core WD 
models by Salaris et al. (2010), He-core tracks by Bedin et 
al. (2008a) and Serenelli et al. (2002), and additional CO-core 
models calculated specifically for this project. We focus here 
on the mF275w versus mF275w - 'Wf438w CMD, which shows a 
very clear separation of the two WDCSs. 

The brightest few magnitudes of a WDCS are populated by 
objects that just entered the cooling sequence and are therefore 
cooling extremely rapidly. We can thus approximate all the 
stars in this transient stage as having the same WD mass, since 
they presumably all started with the same tumoff (TO) mass. It 
is therefore also safe to match them with single-mass cooling 
tracks rather than formal WD isochrones. Theory predicts a 
typical WD mass of ~ O.54A10 in old populations (see, e.g., 
Weiss & Ferguson 2009), which is consistent with what Kalirai 
et al. (2004) found in a spectroscopic study of six bright DA 
WDs in M 4. 

There are a few possible theoretical explanations for the 
split of the bright WDCS of u Cen: (1) unresolved WD-nWD 
binaries, (2) a combination of DA+DB objects, or (3) bimodal 
WD mass distribution, presumably related to the He-enhanced 
subpopulation. 

Regarding the binary scenario, we note that the split WDCS 
at these bright magnitudes cannot be explained by unresolved 
~ O.54AI0 WD+WD binaries. Although the split we observe is 
very nearly equal to that expected for equal-mass binaries, the 
masses of the two components would have to be nearly iden- 
tical at birth, since only an extremely narrow range of initial 
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Fig. 4. — Left panels: one of the observed CMDs of Fig.|3](top) and 
a theoretical counterpart (bottom). WD tracks for CO-core 0.55yV1o 
DA (blue solid) and DB models (blue dash-dotted), 0.46A1o DA CO- 
core (red dotted) and He-core (red dashed) models are fit to the ob- 
served CMD. In the theoretical panels, stars are color coded in azure 
(0.55A1g) and red (0.46 Alo)- Right panels: one of the observed 2CDs 
of Fig. [3] (top) and a theoretical counterpart (bottom). 

mass occupies the transient cooling phase at the same time. If 
the mass ratio of the progenitors is distributed approximately 
homogeneously between and 1 (as found by Milone et al. 
2012), then the two components would be continuously dis- 
tributed over a large magnitude range, and the net effect would 
be a moderate widening of the bright WDCS , not a split (as 
shown for example by the simulations for NGC 679 1 , Bedin et 
al. 2008b). 

We now explore the possibilities that the split is caused by 
either a combination of DA-hDB objects or a WD mass di- 
cothomy. The top-left panel of Fig. |4] superposes on the ob- 
served WDCS four theoretical cooling sequences: a O.SSAI© 
CO-core DA track (blue solid line), the DB counterpart (blue 
dash-dotted line), a 0A6A{q CO-core DA track (red dotted 
line), and a He-core track (and DA atmosphere) with the same 
0.46 Al© mass (red dashed line). 

The O.55AI0 DA model that fits the bWD sequence for the 
eclipsing binary distance represents the "canonical" WD popu- 
lation. The DB track for the same mass is very close to the DA 
track in the region of interest, so the split cannot be explained 
by DA and DB populations with the same mass. We note 
that a bimodal H-layer mass distribution is also ruled out be- 
cause any sequence with lower thickness (the "standard" value 
10 "^AlwD employed in our models is an upper evolutionary 
limit) is shifted to bluer colors and cannot explain the red se- 
quence. 



The only explanation for the rWD objects is therefore a se- 
quence of lower-mass WDs, which would conseguently have 
larger radii. We find that the rWD sequence can be reproduced 
by ~ 0.46 AI0 models, with either CO- or He-cores. 

The bottom left-hand panel of Fig. |4] shows a CMD ob- 
tained from a composite synthetic sample made of the O.SSAl© 
DA and ~ 0.46A1g He-core models. The photometric error 
has been implemented following the results of the data re- 
duction. The synthetic CMD reproduces the bright split se- 
quence, which in the simulation merges around /wpavsw ~ 23.5 
-consistent with the observations- due to the increased photo- 
metric error. The right-hand panels of Fig. H]show one of the 
2CDs of Fig. [3] both observed (top) and simulated (bottom). 
Notice that, like the observations, the simulated sample also 
does not show the bWD+rWD split in the color-color plane. 

5. The Origin of the rWD Sequence 

It is tempting to map the two WD populations to the two 
main MS components, and associate the rWD sequence with 
the helium-rich MS population. To test this hypothesis, we can 
also use the ratio of bWD to rWD stars, and compare it with 
theoretical predictions. First of all, we need to know what kind 
of WDs each is made up of, since CO and He WDs cool at 
different rates. An answer to these questions requires a brief 
examination of the cluster HB morphology. 

Cassisi et al. (2009, C09) showed that the progeny of the 
He-enhanced MS (Y - 0.40) must be confined to the blue tail 
of the observed HB, and found that the blue clump of stars at 
the end of the extended HB tail (see Fig. 2 of C09) can be re- 
produced by either hot flashers or canonical He-rich HB stars, 
or a mixture of both types. 

The minimum mass for the He-rich component along the 
HB is ~ 0.46A4o (C09). These objects will form CO-core WDs 
with essentially the same mass. Even a small amount of addi- 
tional mass loss along the red giant branch (RGB) would pre- 
vent core He-ignition and produce ~ O.46AI0 He-core WDs. 
Note that these stars with high initial He abundances need to 
lose only ~ 0.10-0.15 Mq along the RGB to fall below the 
limit for He-ignition, since their TO mass is about ~ Q.IMq 
less than that of the normal-He stars. 

The presence of the "mass gap" between the two WD se- 
quences (0.55 vs. 0.46 AI0) supports the idea that the blue 
clump in the HB is populated mainly by the evolved mem- 
bers of the He-rich subpopulation, which eventually become 
CO and -possibly- He-core WDs with a mass of ~ O.46AI0, 
creating the rWD sequence. On the other hand, the bWD se- 
quence is instead produced by the end-products of He-normal 
stars (Y - 0.25) that populate the rest of the observed HB. 
Although a fraction of the He-normal population could ex- 
plain the red end of the blue clump, such stars would have 
~ 0.49 AI0 cores and the resulting WDs, with essentially the 
same mass (because of the vanishingly small envelopes of their 
HB progenitors), would occupy the gap between the observed 
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sequences. 

With this understanding of the evolutionary origin of the 
rWD sequence, we investigate whether our interpretation can 
also explain the observed A^rWD/M)WD = 1 .8 + 0.2 number ratio 
in the magnitude range 20.4 < otf275w < 23.4. Our CMDs 
show that the ratio between blue-clump stars and the rest of 
the HB population is ~ 0.32, and the ratio between the He- 
normal and He-rich populations along the MS is ~ 1 in this 
central field. 

We consider two 13.5 Gyr isochrones from the BaSTI 
database (Pietrinferni et al. 2006): one with [Fe/H]=-1.6, 
Y=0.246 to represent the He-normal population, and the other 
with [Fe/H]=-1.3, Y=0.40 to represent the He-rich compo- 
nent. The exact age is not critical: if the age of both compo- 
nents is decreased by 1-2 Gyr and/or the He-rich population is 
made younger by ~ 1 Gyr, the quantitative result of our analy- 
sis is largely unchanged. The isochrones allow us to calculate 
the so-called evolutionary flux, which tells us the number of 
stars leaving the MS per year: 



bit) = O(Mto) 



dt 
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TO 
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(Renzini & Buzzoni 1986). The adopted mass function 
<1>(A1to) has the Salpeter slope. The TO mass A^to and its 
time-derivative are determined from the isochrones, and the 
constant A is fixed by constraining that the number of MS 
stars (below the TO) of the two components to be the same. 
(Note that the MS is the only region that allows us to deter- 
mine the mass-function normalization and initial fraction of 
He-rich stars). The value of b{t) of the He-enhanced popu- 
lation turns out to be 1.1 times the value for the He-normal 
component. Given that beyond the TO (apart from the faint 
end of the WDCS) the initial value of the evolving mass along 
an isochrone is ~ Mjo, the number of stars in a given post-MS 
phase can be approximated by A^pms = b{t) x ^pms (where /pms 
is the lifetime of the post-MS evolutionary stage). 

The typical HB mass for stars in the He-normal component 
(i.e., not in the blue clump) is ~ 0.61-0.62 Mq. Using b{t), 
the model HB timescales provide an expected number ratio 
of ~ 1.7 between blue-clump stars (He-rich) and the rest of 
the HB population. The observed ratio is ~ 0.32, thus only 
~ 20 % of the He-rich population must have produced blue- 
clump stars. This means that ~ 80 % of the He-burning stars 
that we expect from the He-rich component seen along the 
MS, must have missed the He-ignition and evolved onto the 
WDCS, as He-core WDs. If these stars left the RGB close to 
He-ignition (where mass loss is expected to be more efficient), 
they now populate the rWD sequence as ~ O.46AI0 He-core 
WDs. Likewise, the progeny on the (He-rich) HB blue clump 
produced ~ 0.46 A^© CO-core WDs, and the He-normal popu- 
lation produced "canonical" O.SSAIq CO-core WDs. 

We can now split the evolutionary flux for the He-rich pop- 
ulation into two components, according to the fraction of stars 
that achieve He-ignition. Over the magnitude range of interest. 



the 0.46 AI0 CO-core model has approximately the same cool- 
ing time of the 0.55 Al© model (on the order of 10*^ yrs), while 
the 0.46 AIq He-core model has cooling timescales a factor of 
~ 2.4 longer By considering these b{t) and timescale ratios, 
we derive an expected number ratio: 



/7(prog. rWD)"'= '^"'^ x 



A^bWD 



fo(prog. bWD) X fbwD 

£0 
rWD 



/7(prog. rWD)'^° x tf° 



/7(prog.bWD)xfbWD 
and, since fJ?5n°'7fbWD = 2.4 and fS?n"7rbWD = 1-0, we have 



rWD 

A^rWD 
A^bWD 



0.8 X 1.1 X 2.4 + 0.2 X 1.1 X 1.0 ^ 2.3, 



not far from the observed value of 1 .8 + 0.2. This also implies 
that the number of He-core WDs along the rWD sequence is 
^ (0.8 X 1.1 X 2.4)/(0.2 X 1.1 X 1.0) - 10 times the number of 
CO-core WDs. 

To summarize, when taking into account the various evo- 
lutionary paths during the post-MS phases of the He-rich and 
He-normal MS components and the appropriate timescales, we 
can reasonably reproduce the observed A^rWD /A^bWD ratio, and 
derive a very high fraction of He-core WDs along the red CS. 

6. Summary and Conclusions 

Our high-precision, empirical-PSF-based photometric tech- 
niques have revealed a never-before-seen split in the bright part 
of the WDCS in a Galactic GC. We have interpreted the ob- 
served WDCS as follows: 

• The blue sequence corresponds to ~ O.55A10 CO-core 
WDs from the He-normal MS population, and the red 
sequence corresponds to ~ O.46AI0 both CO-core and 
He-core WDs from the He-enriched population. The ra- 
dial gradients of the various populations are consistent 
with this scenario. 

• The "mass gap" between the two WD sequences means 
that the blue clump along the cluster HB has to be pop- 
ulated mainly by stars belonging to the He-rich compo- 
nent. 

• The observed number ratio between blue-clump HB 
stars and the rest of the HB population is consistent with 
the number ratio of He-rich to He-normal MS stars, only 
if ~ 80 % of the HB stars expected to be produced by 
the He-rich component failed to achieve He-ignition and 
evolved onto the red WDCS as ~ 0.46 AIq He-core WDs. 

• It is not extreme mass loss that produces these He-core 
WDs, but rather the fact that the TO mass of the He-rich 
subpopulation is ~ O.lAi^ lower than that of the He- 
normal component. A total mass loss of just ~ O.15A40 
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along the RGB suffices to prevent central He-ignition for 
this subpopulation. 

• The number ratio between CO-core and He-core WDs 
(with longer cooling times) along the red WDCS is 
~ 0.1. This mixture of WDs along the red WDCS 
approximately explains the observed number ratio 
(rWD/bWD). This large contribution of massive He-core 
WDs confirms the more indirect inference by Calamida 
et al. (2008), that the bright WD sequence must be pop- 
ulated by a substantial fraction of He-core WDs. 

In conclusion, our observations and theoretical analysis 
have for the first time disclosed a connection between glob- 
ular clusters with extreme blue HB stars and enhanced He, and 
the morphology - hence mass distribution - of their WDCS. 
It would be extremely important to study the bright WDCS in 
other clusters with extreme HBs and He-rich subpopulations 
(i.e. NGC 6752, NGC 2808) to see whether a spUt WDCS 
analogous to the one detected in w Cen can be found. 
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